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In a recent paper, Luck Q investigated the critical behaviour of one-dimensional Ising 
quantum chains with couplings constants modulated according to general non-periodic 
sequences. In this short note, we take a closer look at the case where the sequences are 
obtained from (two-letter) substitution rules and at the consequences of Luck's results at 
criticality. They imply that only for a certain class of substitution rules the long-distance 
behaviour is still described by the c=l/2 conformal field theory of a free Majorana fermion 
as for the periodic Ising quantum chain, whereas the general case does not lead to a 
conformally invariant scaling limit. 
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1 Introduction 

In this note, we consider the Ising quantum chain with ferromagnetic exchange couphngs 
Ej > which follow an, in general, non-periodic sequence of finitely many different values. 
For convenience, we choose the (constant) transversal field to be equal to one, and consider 
the following Hamiltonians 

of (anti-) periodic approximants (cr^+i = ±crj^) obtained by truncating the sequence of 
exchange couplings to the first A^ elements. Here, a"' denotes the Pauli matrix a" acting on 
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site j of an A^-fold tensor product space. Note that H^"^' commutes with the operator 

N 

Q = li^h Q' = Id, (1.2) 

which has eigenvalues ±1. We denote the projectors onto the corresponding eigenspaces 
(sometimes also called sectors) by 

P± = liId±Q). (1.3) 

In ref. |]l|, the Hamiltonian (|1 . 1| ) is formulated in terms of fermionic operators by means 
of a Jordan- Wigner transformation and diagonalized by a suitable Bogoljubov-Valatin 
transformation. To be more precise, it is in fact the so-called mixed sector Hamiltonian 
H^~^^ defined by 

^(±) _ H^±)p^ + H^^^P^ (1.4) 

which is considered. This is necessary since the Jordan- Wigner transformation of H^'^'^ 
produces non-local boundary terms in the fermionic operators whereas the above Hamilto- 
nians ( |1.4D corresponds to periodic (resp. antiperiodic) boundary conditions in terms of the 
fermions. 

Let us cite three results of ref. ^ which are central to our subsequent discussion. The 
system described by the Hamiltonian (|1.1| ) resp. (|1.4| ) is critical at yU = with 



1 ^ 
^ = j™ TfElogfe) • (1-5) 

The dispersion relation for the low-energy one-particle excitations with energy A takes the 
following form 

A^ = v^q' + f,'), (1.6) 

where q denotes the momentum, and v, the velocity of elementary excitations, is given by 

-1 T N N k 

provided the latter limit exists. 

In what follows, we consider sequences of coupling constants which are obtained by 
substitution rules, focussing on the case of two-letter substitution rules (although most 
properties can be generalized quite easily to the n-letter case). 

2 Substitution Rules and Matrices 

We consider two-letter substitution rules 

e: ? ^ "" (2.1) 

b ^ Wb 



where Wa and Wh are words in a and h (we do not allow for inverses of a or 6 here). If 
one defines multiplication of words through concatenation, the action of q is extended to 
arbitrary words in a and h via the homomorphism property Q{WaWb) = Q{wa)Q{wb)i for 
details see ref. @] and references included therein. To g we associate a 2 x 2 matrix Rg 

' ~ [ MWa) #.K) ) ^^ ^ 

whose elements count the number of a's and 6's in the words Wa and Wb, respectively. Note 
that we use the transpose matrix in comparison with [Q because we need only the statis- 
tical eigenvectors in our discussion. They are then the right-eigenvectors of R^. With this 
convention, one also has Rgoa = Rg ■ Ru- By using the substitution rule g of ( |2.1D itera- 
tively on an initial word 1^(0), say w(0) = a for definiteness, one obtains a sequence of words 
w{n) = g{w{n—l)). Since we are interested in sequences which have a unique limit word w, 
we restrict ourselves to substitution rules where Wa begins with the letter a. In this case, 
the sequence w{n) = g{w{n — 1)) obviously commences with w{n—l) and thus each iteration 
only appends letters to the previous word. The length (i.e., the number of letters) of the 
word w{n) is given by fg{n) defined as follows 

e,{n) = (%^^A = R,e,{n-1), e,(0) = 

f,{n) = e^;\n) + ef^n) . (2.3) 

We denote the eigenvalues of Rg in (|2.2| ) by A^"^^ where A^^-* stands for the Perron- 
Frobenius eigenvalue. The corresponding statistically normalized eigenvector is determined 

by 

The eigenvalue A^^-* determines the asymptotic inflation factor for one substitution, whereas 
Pa {Pb) is the frequency of the letter a (resp. b) in the limit word. Let us now discuss under 
which conditions A*^^-* contains some information about fluctuations. 



3 Fluctuations 

Consider the truncated sequences w^^^ of the first A^ letters of the limit word w obtained 
from the substitution rule (12.11) with initial word a. To measure the fluctuation, we define 




g{N) = #„(w;(^))-p,iV, g^ = g{f,{n)) , (3.1) 

and 

hiN) = max|^(M)|, K = M/.H) • (3-2) 



Note that it does not matter whether we look at fluctuations in the frequency of the let- 
ter a or 6 since #a('W^*-^-') + #fe('U^*'^'*) = N and Pa + Pb = 1- Hence, g{N) just changes 
sign if one replaces a by 6 in Eqs. ( p.l| ) and ( |3.2| ). Of course, if A^+) is not degenerate, 
\imN^^{g{N)/N) = 0. 

The behaviour of g{N) for words of length A^ = /^(n) which correspond to proper (or 
complete) iteration steps is governed by the second largest eigenvalue A*^"^ [Q]. In fact, writing 
the starting vector 6^(0) ( |2.^ ) as a linear combination of the eigenvectors of R^ (|2.2|), one 
easily verifies gn ~ A''"-' . Hence \\^^^\ < 1 implies that gn converges to zero for n -^ oo. On 
the other hand, if lA*-^-*! > 1 then gn in general diverges. In the limiting case of \\^~^ = 1, 
\gn\ is constant and therefore is bounded away from zero and infinity. This observation brings 
along, once again, the concept of Pisot-Vijayaraghavan numbers (PV-numbers for short) |Q. 
They are real algebraic integers -(9 > 1 all algebraic conjugates of which (except i)) lie inside 
the unit circle. If the characteristic polynomial of Rg is irreducible over the integers and if 
the Perron- Frobenius eigenvalue is larger than 1, the PV-property really is what determines 
the conformal nature of the critical point. The same seems still to be true if the characteristic 
polynomial is reducible but all eigenvalues except the largest one lie inside the unit circle 
- in which case we say that the underlying substitution has bounded fluctuation property. 
However, the reducible case requires some care, as we will demonstrate by an example. 

To understand why this is so important, one has to realize that these considerations only 
apply to words which are obtained by proper iteration steps. In between, fluctuations can 
behave quite differently, especially in the so-called marginal case |A(~^| = 1 |l|, |[. There, 
depending on the actual substitution rule (and not on the substitution matrix alone), one 
can have the situation that h{N) is bounded or that h{N) diverges logarithmically with N 
(or, in other words, hn diverges linearly with n). If |A''~-'| > 1, h(N) can diverge like a power 
law. Again, polynomials reducible over the integers are to be treated carefully, in particular 
for generalizations to the n-letter case. 

As an illustrative example, consider the substitution rules that have the substitution 
matrix 

R. = (II) (3.3) 

with eigenvalues A*^+^ = 3 and A^"^ = 1, and Pa = Pb = 1/2- To obtain a unique limit word, 
we want Wa to commence with a, which leaves us with six different substitution rules: 



(3.4) 



Of these, ^5 is special in the sense that it leads to the periodic sequence ababababa . . . which 
means that in this case the fluctuations are certainly bounded (since in any finite part the 
numbers of a's and 6's differ at most by one). As it turns out, this is only true for this special 
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sequence, in all the five other cases /i„ grows linearly with n. More precisely, one observes 



2hr. 



{ n+1 

n 
1 
maxfl, n — 11 



for ^1, Q2, and ^3 
for ^4 
for ^5 
for ^6 



(3.5) 



for n > 0. 

In Fig. 1, we show the different behaviour of the fluctuations g{N) for four typical substi- 
tution rules, namely the Thue-Morse sequence (a), the Silver Mean sequence (b) (sometimes 
also called "Octonacci" sequence), the Period-Doubling sequence (c), and the Binary non- 
Pisot sequence (d). The substitution rules which deflne these sequences together with their 
statistical properties are summarized in Table 1. Fig. 2 shows the quantities \gn\ ( p.l| ) and 
hn ( ^.21 ) for the four different sequences. The marginal case of the Period-Doubling sequence 
(A*-"-* = —1) clearly shows the linear divergence of hn in n whereas \gn\ is constant. 



Table 1: Four typical substitution rules and the statistical properties of the corresponding sequences 



sequence 


substitution rule 


eigenvalues A^''^^ 


values of Pa and pb 
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4 Critical Point and Fermion Velocity 



We now come back to the Hamiltonian ( |1.4] ). We relate the coupling constants Ej to a 
substitution rule g ( p.l|) in the following way: 



(4.1) 



J Ea if j-th letter in w is a 
■^ I £6 if j-th letter in w is b. 

The condition /i = (|1.5|) for criticality now reads 

Pa log{ea) + Pb ^og{eb) = 

which has the following one-parameter solution 0, |^ 

Sa = r-P>' , Eb -- 



r-Pa 



(4.2) 



(4.3) 



where r > is any positive real number. 

As is well known (compare and references therein), the finite-size scaling limit of 
the periodic Ising quantum chain (which corresponds to r = 1) is described by the c= 1/2 
conformal field theory of a free Majorana fermion. To obtain a conformally invariant scaling 
limit one has to have a linear dispersion relation at criticality, which means that the limit 
in Eq. ( |1.7| ) must exist. This in turn is only to be expected if the fiuctuations g{N) remain 
bounded, i.e., hn < S{g) for all n, since at criticality one obtains from Eqs. ( |3.1| ) and ( |4.3| ) 

N 

Y[e, = r-^W (4.4) 



and all products of this type enter in Eq. ( p..7|) . On the other hand, if the limit in Eq. (|L] 
does exist, the scaling limit will be the same as for the periodic quantum Ising chain. 

Although we are not going to use it in this note, we should add that there is a very elegant 
and efficient way to investigate this kind of systems by means of the corresponding trace 
map (compare 0, §] and references therein). In this context, the critical point is obtained 
from a unique one-parameter family of bounded orbits in the accessible phase space region 
of the trace map P . 

5 Fermion and Conformal Spectrum 



We now consider the spectrum of the Hamiltonian ( |1.4D at criticality. It is described in terms 
of A^ fermion frequencies A^ > 0. In Fig. 3, we present the integrated density of the A^ 
(which we devided by the their largest value for convenience) for coupling constants defined 
by r = 2 ( [4.3| ) and several sizes of the chain. In general there are no exact degeneracies in the 
spectrum in contrast to the periodic case. However, as the size of the system increases, the 
frequencies tend to accumulate which creates the nearly vertical steps in Fig. 3, especially 
close to the maximal frequency. The plots show characteristic gaps in the fermion spectrum 
the locations of which (on vertical axis) are in accordance with the general gap labeling 



theorem [|l^, |TI|, |J]. We should mention though that the Thue- Morse chain does not show 



closed gaps here, in contrast to the situation with electronic spectra of Schrodinger operators 



r^j , a phenomenon that deserves further exploration. 

For models whose continuum limit is described by a conformal field theory, conformal 
invariance specifies the behaviour of the low-energy excitations in the infinite size limit 
N —>■ oo. Essentially, they have to show a leading 1/A^ behaviour and the level and degener- 
acy structure of the spectrum (after appropriate overall scaling of the gaps) is described by 
representations of the Virasoro algebra with central extension c which is the central charge 
of the conformal field theory (see e.g. |]^, [1^, 0). This of course means that it is the lower 
part of the integrated fermion density shown in Fig. 3 which is important for the conformal 
spectra. 



To have a closer look at this part, we consider the scaled energy gaps 

E] = ^iE,-Eo) (5.1) 

of the Hamiltonian H^^^ ( |1.4|) , where Eq denotes the ground-state energy and v is the velocity 
of elementary excitations. For the Thue-Morse and the Silver Mean sequences, v is finite (in 
the limit N -^ oo and for finite r > 0) and given by [|l| 

- — — ' for the Thue-Morse chain 



' ' (5.2) 



log(r) for the Silver Mean chain. 

In fact, the second value generally applies for all quasiperiodic sequences Jll which can 
be obtained from a certain section through a higher-dimensional periodic structure via the 



dualization method [|T5|. On the other hand, the symmetric random dimer chain yields 
the same result as the Thue-Morse chain 0. It thus represents a disordered model (with 
bounded fluctuations) which nevertheless shows the same critical behaviour as the periodic 
Ising chain and, in particular, leads to a conformally invariant continuum limit. 



From our discussion of the fluctuations and from Eq. ( |1.7D above it is clear that for 
the Period-Doubling and Binary non-Pisot sequences the fermion velocity v should vanish. 
In Fig. 4, the scaled spectrum of low-energy excitations is shown for our four exemplary 
sequences. For the Thue-Morse and Silver Mean chains, the fermion velocity (|5.2|) is used 
whereas for the other two sequences we simply normalize the first gap to 1/2. (We did not 
systematically investigate the scaling laws here, although this might give independent access 
to the critical exponents calculated in ||l|). The formation of the so-called conformal towers 
can be seen clearly for the spectra of the Thue-Morse and Silver Mean chain (a look at the 
actual data confirms that the degeneracies are those predicted from conformal invariance), 
whereas the spectra of the other two systems do not show any apparent regularities. As 
expected, the normalization factors grow with the size of the system in these cases. As a 
consequence, the nature of the critical point is quite different, compare [Q], and conformal 
invariance is lost. 

6 Concluding Remarks 

The interest in conformally invariant phase transitions and the almost immediate study of 
quasiperiodic systems of the Fibonacci type (see e.g. 0, |T6|, |l^, |l^) originally has led to the 



somewhat misleading conclusion that conformal invariance is robust with respect to many 
sorts of order or even disorder. Although this is true for quasiperiodic Ising quantum chains 
(because that implies the PV-property or the bounded fluctuation property [0), many other 



ordered structures can be defined which destroy conformal invariance at the critical point. 
Even more, one can consider the bounded fiuctuation situation to be somewhat exceptional 



(e.g., PV-numbers are nowhere dense in [1, cxd), see ||2^) and thus conclude that the general 



case of a non-periodic Ising quantum chain does not lead to a conformally invariant scaling 
limit. 
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Figure 1: Fluctuations g{N), compare ( |3.l[ ), for the following sequences: (a) Thue-Morse, (b) Silver 
Mean, (c) Period-Doubling, and (d) Binary non-Pisot. The lengths of sequences which 
correspond to complete iteration steps are indicated. 
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Figure 2: Fluctuations |g„| (3J^) (Figure 2A) and /i„ ( |3.2| ) (Figure 2B) for the following sequences: 
(a) Thue-Morse, (b) Silver Mean, (c) Period-Doubling, and (d) Binary non-Pisot. 
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Figure 3: Integrated density of normalized fermion frequencies for the Hamiltonian H^^' ( |1.4| ). The 
coupling constants are given by Eq. (4.3) with r = 2 and (al)-(a3) the Thue- Morse sequence 
with n = 5,6,7, (bl)-(b3) the Silver Mean sequence with n = 4,5,6, (cl)-(c3) the 
Period-Doubhng sequence with n = 5, 6, 7, and (dl)-(d3) the Binary Non-Pisot sequence 
with n — 4:,5, 6, where n denotes the number of iterations. The spectra do not have true 
degeneracies, and the integrated density reaches 1 in all cases shown. 
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Figure 4: Scaled low-energy spectra of iJ^^^^ ( |l.4| ) for couplings obtained from Eq. (K^) with r = 2 and 
(al)-(a3) the Thue-Morse sequence with n = 6,7,8, (bl)~(b3) the Silver Mean sequence 
with n = 5,6,7, (cl)-(c3) the Period-Doubhng sequence with n = 6,7,8, and (dl)-(d3) 
the Binary Non-Pisot sequence with n = 5,6, 7, where n denotes the number of iterations. 



Normalization of (a) and (b) is taken from (5.2), while for (c) and (d) the first gap is 
normalized to 1/2. 
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